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Abstract 



After a brief review of galactic cosmic rays in the GeV to TeV energy range, we 
describe some current problems of interest for particles of very high energy. Par- 
ticularly interesting are two features of the spectrum, the knee above 10 15 eV and 
the ankle above 10 18 eV. An important question is whether the highest energy par- 
ticles are of extra-galactic origin and, if so, at what energy the transition occurs. 
| A theme common to all energy ranges is use of nuclear abundances as a tool for 

CO ■ understanding the origin of the cosmic radiation. 

O 



1 Introduction 



The cosmic-ray spectrum falls steeply, decreasing by approximately a factor of 
50 per decade increase in energy when plotted as E<p(E) = dN/dln(E). In the 
^ ■ lowest energy region the flux is high enough so that the elemental composition 

of the primary cosmic-ray nuclei can be studied by direct observations with 
detectors lifted above the atmosphere by balloons or spacecraft. In the GeV 
range, even individual isotopes can be resolved. In this lower energy range we 
have the most detailed information on which to base a model of the origin of 
cosmic rays. 

Above 10 14 eV, where the flux falls below several particles per square meter 
per day, direct measurements of the primaries are no longer practical. On the 
other hand, this energy is large enough so that secondary cascades penetrate 
with a footprint large enough to be measured by an array of detectors on the 
ground. Such an extensive air shower (EAS) array typically has dimensions 
of a fraction of a square kilometer or more and can be operated for years 
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Fig. 1. Summary of the cosmic-ray spectrum. The five lines below 100 GeV indicate 
the individual contributions of five individual groups of nuclei (protons, helium, 
CNO, Mg-Si and Fe). For clarity only subsets [1,2,3,4] of the data are shown. 

rather than days or weeks. An EAS experiment in effect uses the atmosphere 
as a calorimeter, so the showers are classified by total energy per particle 
rather than by energy per nucleon as at low energy. Because of the sparse 
quality of the information, however, the energy per shower is only determined 
with relatively large uncertainty, and the best one can do is to determine the 
relative contributions of groups of elements. The situation is exacerbated by 
large shower-to-shower fluctuations coupled with less than full knowledge of 
the properties of the hadronic interactions that govern shower development. 
The latter problem becomes more severe in the highest energy range, which 
is beyond the reach of current accelerators. 

Figure 1 is a schematic plot of the cosmic-ray spectrum as a function of total 
energy per nucleus. The plot is divided into three regions by the shading. The 
relatively detailed measurements of spectra of individual elements (or groups 
of elements) below 10 14 eV supports a fairly well-developed model of the origin 
of these particles. In particular, the observed energy density and propagation 
time of the cosmic-rays together determine the power required by the sources. 
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This connection suggests an approach that may be helpful in understanding 
the origin of the higher energy particles. The second region (dark shading) 
includes the knee where there is a steepening and perhaps some structure in 
the spectrum that needs explanation. The highest energy region includes the 
ankle, a flattening of the spectrum above 10 18 eV that may be related to a 
transition from particles of galactic origin to those accelerated in extra-galactic 
sources. The question of whether there is a transition to cosmic rays of extra- 
galactic origin, and if so where it occurs, is of considerable interest. Finally, 
the question of whether the spectrum extends beyond 10 20 eV is currently the 
foremost problem in high-energy particle astrophysics because of the difficulty 
of accounting for the absence of energy loss by protons in the microwave 
background radiation if the spectrum is not suppressed above this energy. 

In what follows we discuss the three energy regions in order. A common theme 
is use of the primary composition as a clue to the nature of the sources of the 
particles. 



2 Galactic cosmic rays E < 10 eV 



A classic problem in nuclear astrophysics is the determination of the cosmic- 
ray source abundances [5,6]. This requires a model of the galaxy, including the 
spatial distribution of cosmic-ray sources; density, composition and ionization 
state of the interstellar medium; and the strength and topology of the mag- 
netic field. Important recent investigations of cosmic-ray propagation from 
two different viewpoints may be found in Refs. [7] and [8]. Here we give only 
the main points in a simplified form. The starting point is a determination 
of the abundances of different nuclei near the solar system (but outside the 
heliosphere) . Figure 2 shows the elemental abundances in the cosmic radia- 
tion. The overabundance by several orders of magnitude of secondary elements 
such as lithium, beryllium and boron relative to their general abundance in 
solar system material is a consequence of spallation of the more abundant 
primary nuclei, especially carbon and oxygen. A similar situation occurs for 
the secondary nuclei just below iron. 

Given a knowledge of the spallation cross sections, one has to solve a set 
of coupled propagation equations for the diffusion of the cosmic-rays in the 
turbulent, magnetized interstellar medium. A simplified version of the relevant 
set of equations is 

A:[t ' = Qi(E) - (pcn H a t + —) N t (E) + (3cn H J2^i N k (E). (1) 
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Nuclear abundance: cosmic rays compared to solar system 
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Fig. 2. Comparison of Solar system [9] and cosmic-ray elemental abundances. Nu- 
clear abundances are from [10]; protons and helium are from [11,12]. 

Here Ni(E) is the spatial density of cosmic-ray nuclei of mass i, and % is the 
number density of target nuclei (mostly hydrogen) in the interstellar medium, 
Qi(E) is the number of primary nuclei of type i accelerated per cm 3 per 
second, and cr* and a^-H are respectively the total and partial cross sections 
for interactions of cosmic-ray nuclei with the gas in the interstellar medium. 
The second term on the r.h.s. of Eq. 1 represents losses due to interactions 
with cross section cr, and decay for unstable nuclei with lifetime Tj. 



Essentially, in the simplest diffusion model, the secondary to primary ratios 
determine the product of density of the interstellar medium times the charac- 
teristic propagation time of cosmic-rays before they escape from the Galaxy 
(r esc ). This can be seen by solving Eq. 1 for a secondary nucleus S while 
neglecting its losses during propagation and assuming that Qs ~ 0. Repre- 
senting the parent primary nuclei by a single contribution Np then leads to a 
simplified equation 



Nc 



p/3cr e 



5 g/cm 2 . 



(2) 
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The approximate value given in Eq. 2 requires accounting for several terms 
with appropriate spallation cross sections for its derivation from the data. 
The numerical value in Eq. 2 would imply T esc ~ 3 x 10 years if the average 
density p corresponds to one hydrogen atom per cubic centimeter as in the 
disk of the galaxy. 

Two details of the data are particularly important. One is the ratio Rio of 
the unstable isotope 10 Be to stable 9 Be. 10 Be is unstable to /3-decay with a 
half-life of r 10 = 1.5 x 10 6 years. The two isotopes of beryllium are produced 
in comparable amounts by spallation of heavier nuclei. If Rio as measured is 
comparable to the production ratio, then r esc <C t w ; if Rio ~ then r esc ^> r w . 
The data indicate 

r esc f=s 2 x 10 7 years. (3) 

Given the mean density of the interstellar medium, the relatively large value 
of r esc implies that cosmic ray nuclei spend significant time diffusing in low- 
density regions, possibly above the gaseous disk in the galactic halo, before 
escaping into inter-galactic space. 

The other significant feature is that the ratio of secondary to primary nuclei de- 
creases with increasing energy. From Eq. 2 this implies that r esc decreases with 
energy. This behavior is attributed to energy-dependent diffusion, whereby 
higher energy particles diffuse out of the galaxy more quickly than those of 
lower energy. It should apply to all cosmic rays, including primary protons, as 
well as secondary nuclei if the model is consistent and complete. The observed 
primary spectrum (f>(E) and the source spectrum Q(E) are connected by a 
relation of the form 



with 5 ~ 0.6. Since the observed differential energy spectrum is proportional 
to E~ 2 - 7 , the inferred source spectrum would be oc E~ 21 . 

An inferred spectral index at the source of -2.1 is often cited as evidence 
for first order diffusive shock acceleration, which, in the test-particle approx- 
imation predicts a value of -2.0. The situation is actually more complicated. 
On the one hand, the relation 5 cannot hold over a large energy range with- 
out coming in conflict with the observed isotropy of high-energy cosmic-rays. 



c/>(E) = Q(E) x r esc (E). 



(4) 



A simple power-law fit to available data gives 



r esc (E) oc E- s 



(5) 
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Extrapolating Eq. 5 to 10 15 eV, for example, would lead to a value of r esc al- 
most as small as the light travel time across the galactic disk, implying a larger 
anisotropy than is observed. On the other hand, a more complete treatment of 
shock acceleration that accounts for the non-linear back-reaction of the accel- 
erated particles on the shock structure [13], gives a somewhat concave energy 
spectrum rather than a single power law (i.e. steeper spectrum at low energy, 
hardening at high energy). Models in which some of the energy-dependence of 
the secondary /primary ratio is attributed to reacceleration [14] during prop- 
agation can accommodate a somewhat steeper average source spectrum and 
correspondingly a larger value of r esc at high energy, more consistent with 
the observed small anisotropy. One possibility is that the observed average 
spectrum is a composite of contributions from many individual sources having 
a range of individual properties. In what follows, for illustration we assume 
Q{E) oc E~ 2A in which case 5 « 0.3. 

The next step is to estimate the power required of the sources to maintain 
the observed spectrum in equilibrium. The total source power is related to the 
locally observed energy spectrum by 



The pre-factor on the r.h.s. of this equation converts the measured cosmic- 
ray flux (j) to the corresponding spatial density of cosmic rays. Fixing r esc in 
the GeV range according to Eq. 3 and evaluating the integral using the mea- 
sured spectrum and Eq. 5 for r esc leads to an estimate of the source power 
as ~ 10 -26 erg cm _3 s _1 . If the flux observed locally near the solar system is 
typical of the galactic disc, then the total power of the cosmic-ray sources is 
obtained by multiplying by the volume of the disk, ~ 10 67 cm 3 , for an estimate 
of 10 41 erg/s. The average power generated in kinetic energy of ejected mate- 
rial (exclusive of neutrinos) in galactic supernova explosions is ~ 10 42 erg/s 
(10 51 erg/30 year). 

The combination of an acceleration theory which produces approximately the 
right spectrum with high efficiency (first order diffusive shock acceleration) 
and a source with the right power makes supernova shock acceleration a fa- 
vored candidate for the origin of cosmic rays in the galaxy. From observations 
of synchrotron radiation, it is possible to see direct evidence of electrons be- 
ing accelerated to high energy in individual supernova remnants [15]. What 
has proved more difficult is to find examples of individual supernova rem- 
nants in which acceleration of protons can be verified by the production of 
neutral pions [16]. The problem is that it is usually possible to explain the 
observed electromagnetic radiation as radiation from energetic electrons, by 
bremsstrahlung at lower energy and by inverse Compton scattering at higher 
energy. The characteristic kinematic signature of pion production (a peak at 




(6) 
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a photon energy of half the pion mass) tends to be obscured by superimposed 
synchrotron and bremsstrahlung photons [17]. Indirect arguments related to 
the shape of the observed gamma-ray spectrum above a TeV can in some cases 
most easily be accounted for if the source is production of 7r° by protons rather 
than radiation by electrons [18]. A nice review of models of galactic cosmic 
ray sources is Ref. [19]. 

Observation of neutrinos from SNR (or any other potential cosmic accelera- 
tor) would be conclusive evidence for acceleration of primary protons because 
only hadronic processes can produce neutrinos. The expected fluxes are low, 
however, so large detectors will be needed to detect them [20]. 



3 The knee: 10 14 < E < 10 17 eV 

Diffusive shock acceleration works to the extent that charged particles gain 
energy by an amount AE proportional to E each time they cross from up- 
stream to downstream and back upstream of the shock. (Here upstream is the 
unshocked region and E is the energy before adding the increment AE.) After 
a time T the maximum energy achieved is 

E max ~ Zej3 s x B xTV s , (7) 

where (3 S = V s /c refers to the velocity of the shock. This result is an upper 
limit in that it assumes a minimal diffusion length equal to the gyroradius of 
a particle of charge Ze in the magnetic field B upstream of the shock. Using 
numbers typical of Type II supernovae exploding in the average interstellar 
medium gives E max ~ ZxlOO TeV [21]. More recent estimates give a maximum 
energy larger by as much as order of magnitude or more for some types of 
supernovae [22]. 

The nuclear charge, Z, appears in Eq. 7 because acceleration depends on 
the interaction of the particles being accelerated with the moving magnetic 
fields. Particles with the same gyroradius behave in the same way. Thus the 
appropriate variable to characterize acceleration is magnetic rigidity, R = 
pc/Ze ~ E tot /Ze, where p is the total momentum of the particle. Diffusive 
propagation also depends on magnetic fields and hence on rigidity. For both 
acceleration and propagation, therefore, if there is a feature characterized by 
a critical rigidity, R*, then the corresponding critical energy per particle is 
E* = Z x R*. 

The knee of the spectrum is the steepening that occurs above 10 15 eV, as shown 
in Fig. 3, while the ankle is the hardening around 3 x 10 18 eV. One possibility 
is that the knee is associated with the upper limit of acceleration by galactic 
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Fig. 3. High-energy cosmic-ray spectrum. References to the data are given in [23]. 
The shaded region indicates a factor a 30 in total energy (see text). 

supernovae, while the ankle is associated with the onset of an extragalactic 
population that is less intense but has a harder spectrum that dominates at 
sufficiently high energy. A generic two-component model was suggested long 
ago by Peters [24], but this explanation does not fit the data, at least not in 
its simplest form. 

Consider first the situation if all galactic sources accelerated particles to the 
same maximum rigidity, say 10 15 V. Then, since the abundant nuclei are in 
the range 1 < Z < 26, the spectrum would be expected to cutoff altogether 
within a factor of about 30 as a consequence of the prefactor Z in Eq. 7. 
The protons would cut off first, followed by helium then CNO, etc. [24]. The 
range over which the cutoff would occur is indicated by the shaded region 
in Fig. 3. Clearly, this is not at all what happens. Instead, the spectrum 
continues smoothly for another two decades in energy. Even postulating a 
significant contribution from elements heavier than iron (up to uranium [25]) 
cannot explain the smooth continuation all the way up to the ankle. 

One possibility is that most galactic accelerators cut off around a rigidity of 
perhaps 10 15 eV, but a few accelerate particles to much higher energy and 
account for the region between the knee and the ankle [26]. This scenario 
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would be a generalization of Peters' model. Its signature would be a sequence 
of composition cycles alternating between light and heavy dominance as the 
different components from each source cut off. As emphasized by Axford [27], 
however, the problem with this type of model is that it requires a fine-tuning 
of the high-energy spectra so that they rise to join smoothly at the knee then 
steepen to fit the data to ~ 10 18 eV. As a consequence, several models have 
been proposed in which the lower-energy accelerators (E < 10 15 eV) inject 
seed particles into another process that accelerates them to higher energy. 
In this way the spectrum above the knee is naturally continuous with the 
lower energy region. Groups of supernovae [27] and a termination shock in the 
galactic wind [28] have been suggested [29]. 

The fine tuning problem (i.e. to achieve a smooth spectrum with a sequence 
of sources with different maxima) was actually clearly recognized by Peters 
in his original statement of this idea [24]. He correctly point out, however, 
that since the cutoff is a function of rigidity while the events are classified by 
a quantity close to total energy, the underlying discontinuities are smoothed 
out to some extent. An interesting question to ask in this context is what 
source power would be required to fill in the spectrum from the knee to the 
ankle. The answer depends on what is assumed for the spectrum of the sources 
and the energy dependence of propagation in this energy region. Reasonable 
assumptions (e.g. Q(E) oc E~ 2 and r esc oc E~ s with 5 ~ 0.3) lead to an 
estimate of ~ 2 x 10 39 erg/s, less than 10% of the total power requirement 
for all galactic cosmic-rays. For comparison, the micro-quasar SS433 at 3 kpc 
distance has a jet power estimated as 10 39 erg/s [30]. 

Another possibility is that the steepening of the spectrum at the knee is a 
result of a change in properties of diffusion in the interstellar medium such 
that above a certain critical rigidity the characteristic propagation time r esc 
decreases more rapidly with energy. If the underlying acceleration process 
were featureless, then the relative composition as a function of total energy 
per particle would change smoothly, with the proton spectrum steepening 
first by 0.3, followed by successively heavier nuclei. It is interesting that this 
possibility was also explicitly recognized by Peters [24]. 

A good understanding of the composition would go a long way toward clarify- 
ing what is going on in the knee region and beyond. A recent summary of direct 
measurements of various nuclei shows no sign of a rigidity-dependent compo- 
sition change up to the highest energies accessible (~ 10 14 eV/nucleus) [31]. 
The change associated with the knee is in the air shower regime. Because of 
the indirect nature of EAS measurements, however, the composition is diffi- 
cult to determine unambiguously. The composition has to be determined from 
measurements of ratios of different components of air showers at the ground. 
For example, a heavy nucleus like iron generates a shower with a higher ratio 
of muons to electrons than a proton shower of the same energy. The best indi- 
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cation at present comes from the Kascade experiment [32] , which shows clear 
evidence for a "Peters cycle" , the systematic steepening first of hydrogen, then 
of helium, then CNO and finally the iron group. The transition occurs over an 
energy range from approximately 10 15 eV to 3 x 10 16 eV, as expected, but the 
experiment runs out of statistics by 10 17 eV, so the data do not yet discrimi- 
nate among the various possibilities for explaining the spectrum between the 
knee and the ankle. 



4 The ankle: E > 10 17 eV 



The energy range above 10 is where a transition from cosmic rays accelerated 
in the Galaxy to cosmic rays from extra-galactic sources may be expected. The 
gyroradius of 10 19 eV protons is of order 10 kpc, exceeding the dimensions of 
the inner Galaxy. Such particles could not be contained in the Galaxy [33], 
not to mention the difficulty of accelerating particles to such high energy with 
known potential galactic accelerators. The spectrum of the highest energy 
cosmic rays indeed shows an 'ankle' around 3 x 10 18 eV where the power law 
spectrum flattens to a differential index of about 2.7. One possibility is that 
this feature itself reflects the transition to dominance of extragalactic cosmic 
rays. Another possibility is that particles of extragalactic origin become the 
dominant population at lower energy and that the ankle is a feature due to 
energy loss of protons to electron-positron pair production on the microwave 
background during propagation [34]. 

What is needed to decide the question is knowledge of the composition as 
a function of energy. The idea is that in the range where the last galactic 
accelerator reaches its upper limit the composition will start changing in the 
same way as it does at the 'knee'. In an energy range of about one and a half 
decades the composition should change from one dominated by protons and 
light nuclei to heavy nuclei, signaling the upper limit of Galactic sources. As 
the extragalactic population of particles begins to dominate, the composition 
should change back toward lighter nuclei, assuming the extragalactic cosmic 
rays consist primarily of protons. 

In this energy range the composition is measured by the energy dependence 
of the position of shower maximum, X max . An air shower consists of a su- 
perposition of electromagnetic cascades initiated by photons from decay of 7r° 
particles produced by hadronic interactions along the core of the shower as 
it passes through the atmosphere. Most of the energy of the shower is dissi- 
pated by ionization losses of the low-energy electrons and positrons in these 
subshowers. The composite shower reaches a maximum number of particles 
(typically 0.7 particles per GeV of primary energy) and then decreases as the 
individual photons fall below the critical energy for pair production. Because 
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Fig. 4. Plot of data on mean depth of maximum vs energy. Filled squares are data 
from the stereo Fly's Eye [39]. Open symbols show the data of HiRes (squares) and 
HiRes prototype with MIA (circles) [40]. 

each nucleus of mass A and total energy Eq essentially generates A subshowers 
each of energy Eq/A the depth of maximum depends on E /A. Since cascade 
penetration increases logarithmically with energy, 



where Xer is a parameter (the "elongation rate" ) that depends on the under- 
lying properties of hadronic interactions in the cascade. 

All contemporary interaction models predict an elongation rate between 50 
and 60 g/cm 2 per decade of energy for protons. Showers generated by iron 
nuclei would have a similar elongation rate, but would have X max shallower 
by 90 to 100 g/cm 2 at the same total energy (see Eq. 8). In general, the depth 
of maximum should reflect a changing composition according to Eq. 8. An 
extreme case would be a composition changing from pure iron at the end of 
the galactic component to pure proton at higher energy if all extragalactic 
cosmic rays are protons. In this case X max would increase by « 100 g/cm 2 
more than the amount due to the elongation rate alone in the energy interval 
in which the transition occurs. 

Atmospheric fluorescence telescopes such as the Fly's Eye [35] directly mea- 



X, 



max 
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sure the longitudinal development of large air showers. After correcting for 
atmospheric absorption and Cherenkov light and accounting for energy lost 
to neutrinos and energetic muons, the measured shower profile can be fit- 
ted and a value of X max determined for each shower [36]. Figure 4 shows 
the average depth of shower maximum as measured by three fluorescence ex- 
periments. The lines show X max predicted by different hadronic interaction 
models [37,38]. The model dependence is not high for iron initiated showers, 
but reaches about 30 g/cm 2 for very high energy proton showers. 

Taken together, the data show a transition from a large fraction of heavy nuclei 
around 10 17 eV toward the proton predictions above 10 18 5 eV. Whereas the 
original Fly's Eye stereo measurements [39] suggest a mild transition toward 
protons between 10 17 ' 5 and 10 19 eV, the more recent HiRes/MIA prototype 
data [40] data suggest that the transition begins already at 10 17 eV and is 
complete by 10 18 eV. The uncertainty about the composition in this region is 
indicative of the difficulty of the measurements compounded by uncertainties 
in the hadronic interaction models used to interpret the data. Alan Watson 
has given a nice review of these problems recently [41]. He points to progress 
in solving them as a key to understanding the origin of the ultra-high energy 
cosmic radiation. 



25.5 
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Fig. 5. Energy spectrum of the cosmic rays detected by the Fly's Eye (filled cir- 
cles) [35]. For extension to lower energy the spectrum measured by the Akeno de- 
tector is shown with filled squares [2]. 

Finally we turn to the energy spectrum itself in Fig. 5. Rather than showing all 
data sets as in Fig. 3, we show only the Fly's Eye stereo data at high energy. 
The stereo data consist of air showers detected simultaneously by the two 
Fly's Eye fluorescence detectors. Such events are better reconstructed and have 
smaller systematic uncertainties than showers detected by a single fluorescence 
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telescope. However, the exposure for stereo observations is significantly smaller 
and statistical uncertainties become very large above 10 19 eV. The advantage 
of looking at a single data set is that spectral features become clear. The 
spectrum appears to steepen slightly for E > 10 17 ' 5 eV, a feature referred to 
as the "second knee". This is followed by the flattening around 10 18 5 eV which 
is the ankle. Spectra measured by other experiments [3,4,42] show qualitatively 
similar features, except for an apparent excess of events with E > 10 20 eV in 
the AGASA experiment [3]. 

If the second knee is associated with the end of the galactic cosmic rays, then 
the ankle would most naturally be attributed to pair losses during propagation 
in the microwave background, as proposed by Berezinsky [34]. In this case, the 
cosmic rays down to ~ 10 17 ' 5 eV could be primarily of extragalactic origin. If 
instead the ankle is where the extra-galactic component becomes the dominant 
population, then most of the particles could be of galactic origin up to ~ 
10 18 - 5 eV. 

The importance of this difference becomes clear when we take the next step 
and estimate the power needed to supply the extra-galactic component. For 
a cosmological distribution of sources, the analog of r esc in Eq. 1 is the age 
of the Universe, r# ~ 10 10 yrs. The required power is then given by Eq. 6 
with r esc — > Tfj and with the total flux replaced by the extragalactic com- 
ponent, 4>eg- The result therefore depends on where the extra-galactic com- 
ponent is normalized to the observed spectrum and on how it is extrapo- 
lated to low energy below the galactic component. Assuming a differential 
spectral index of —2 and normalizing at 10 19 eV, the estimated energy den- 
sity is 2 x 1CT 19 erg/cm 3 . This leads to an estimated power requirement of 
~ 10 37 erg/Mpc 3 /s. Both gamma-ray burst sources (GRB) [43,44] and active 
galactic nuclei (AGN) [34] are leading contenders for sources of the ultra-high 
energy cosmic rays, in part because their estimated intrinsic luminosities sug- 
gest that they are powerful enough to satisfy this requirement, given their 
observed distributions. Shifting the normalization point down in energy to 

10 18 eV requires a factor 10 more power, which is probably easier to satisfy 
with AGN. In the GRB model the extragalactic component is normalized at 

10 19 eV [45]. 



5 Is there a GZK cutoff? 



This may be the most discussed question in particle astrophysics. It refers to 
the fact that, for a cosmological distribution of sources, protons with energy 
above threshold for photo-pion production on the microwave background are 
expected to lose energy to this process [46,47]. Nuclei also lose energy by 
photo-disintegration. The result should be a suppression of the flux for energies 
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above 5 x 10 19 eV. The two experiments with the largest current exposure [3,4] 
show results which give a different indication about whether the expected 
suppression is there or not. This can be seen at the high end of Fig. 1. The 
main problem is the extremely low flux [48], which is of order one particle 
per square kilometer per century above 10 20 eV, coupled with the limited 
exposures of the detectors. Answering this question is a major focus of the 
Auger Project [49]. 

If the spectrum continues past the GZK energy without suppression, high- 
statistics studies will be needed to look for clustering or anisotropies related 
to individual sources, which would presumably be nearby in the absence of 
a cutoff. If a suppression is observed, then accumulation of statistics would 
be useful to understand the shape of the spectrum around the GZK feature, 
which carries information about the distribution and evolution of sources on 
cosmological scales [50]. If the spectrum of an extra-galactic component can be 
sufficiently well determined also at lower energy by understanding the sources, 
then it could be subtracted from the total observed spectrum revealing the end 
of the spectrum of galactic cosmic rays. Preliminary and speculative versions 
of this exercise appear in Ref. [34] and Ref. [51]. 
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